o 



in 

o 



0^ 



X 



Isospin Dependence of Proton and Neutron Radii 
within Relativistic Mean Field Theory * 

M. Warda, B. Nerlo-Pomorska and K. Pomorski 
Theoretical Physics Department, University M.C.S., Lublin, Poland 

Abstract 



The binding energies, shapes and sizes of even-even /3-stable nuclei with ^ > 40 
^ ' and a few chains of isotopes with Z=50, 56, 82, 94 protons and isotones with A^=50, 

^«0 . 82, 126 neutrons are analyzed. The average isospin dependence of the radii of pro- 

tons and neutrons evaluated within the relativistic mean field theory is studied. A 



Psj ' simple, phenomenological formula for neutron radii is proposed. 
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C ■ 1 Introduction 



The nuclear radii are ones of the crucial quantities good for testing of every theoretical 
model of nucleus. The radii are measured with a high accuracy for the charge density 
distributions [^ and less precise for the neutron ones 0. The data concerning nuclear 
sizes and shapes are already known for many nuclei, especially for those close to the /3 
stability line, but also more exotic nuclei are explored extensively by experimentalists at 
present. It would be worthwhile to know what the relativistic mean field theory (RMFT) 
predicts for such nuclei and to approximate the results in terms of a simple, practical 
formulae for radii. 

The nuclear radii depend chiefly on the number of nucleons {A) 

Ro = roAV3 , (1) 

where the radius constant tq ~1.2 fm. It is the result of saturation property of nuclear 
forces which is manifested in the experimental fact that the volume of nucleus is roughly 
proportional to the mass number A. However, the formula (1) is not valid any more for 
nuclei in which numbers of protons (Z) and neutrons (A^) differ significantly. 

*This work is partly supported by the Polish Committee of Scientific Research under Contract No. 
2 P03B 049 09. 



It was found in Ref. [Q that the isospin dependent formula for the nuclear charge 
radius constant 

rf = 1.25 (l - 0.2 ■ ^^^) fm, (2) 

describes much better than Eq. (1) the experimentally known charge mean square radii 
of even-even nuclei with A > 60. The formula (2) was obtained assuming the uniform 
charge distribution within deformed nucleus. In such an approximation the root mean 
square radius (RMSR) of deformed nucleus is given by the following formula 

<r'>'/'=^Ro-g{e,e,), (3) 

where the function g describes the dependence of the mean square radius on deforma- 
tion. The equilibrium deformations of nuclei were taken from Ref. [^, where the two- 
dimensional {e, £4) space of deformation parameters was used. The potential energy 
surfaces were calculated in Ref. by the Strutinsky prescription with the zero-point 
energy correction terms according to the generator coordinate method f^. 

Lateron it was found, after more broad calculations with the equilibrium deforma- 
tions taken from Ref. 0, that the additional term ~ ^ in the formula (2) should appear 
when reproducing the mean square radii of all nuclei beginning from the lightest ones up 
to the actinides [§] 

r[)'' = 1.240 (l - 0.191 -— + 1.646 ■ -\ fm . (4) 

A further development of the formula for nuclear radius was made in Ref. [Q after 
extensive Hartree-Bogolubov calculations with various sets of Skyrme effective interactions 
for proton and neutron density distributions in a few spherical nuclei. The authors of Ref. 
postulate new terms proportional to l/A"^ in the formula for Tq^. 

The aim of our present work is to find simple formulae which approximate the results 
obtained for proton and neutron radii within the RMFT pD| , |TT|. We investigate all even- 
even P stable nuclei as well as the isotopes and isotones corresponding to major magic 
Z and N numbers. When discussing the nuclides far from stability we change the mass 
number A up to the zero neutron or proton separation energy. 

In Section 2 we describe briefly the theoretical model and its parameters. In Section 3 
we analyze the results of the calculations and we write formulae for the radii of the 
proton and neutron distributions, which approximate the theoretical results of the RMFT. 
A simple, phenomenological formula for neutron radii is proposed as the result of the 
above analysis. The conclusions and further investigations perspectives are described in 
Section 4. 



2 Theoretical model 



The relativistic mean field theory |[T^ is a variational model based on a standard La- 
grangian density []TT| 
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consisting of nucleon tp, mesons cr,uj,p and electromagnetic A fields. The a mesons 
potential has been taken in the nonlinear form: 
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(6) 



It was found in |]T2| that the NL-3 |T3[ parameters set of the mean field Lagrangian 
(5) reproduced well binding energies, proton and neutron separation energies, electric 
quadrupole moments and radii of all nuclei along the whole /3-stability line. The NL-3 
parameters are: 

- nucleon mass M = 939 MeV 

- meson masses m„ = 508.194 MeV, m,, = 782.501 MeV, 



rrir 



763 MeV 



- meson coupling constants g^ = 10.217, g^= 12.868, gp = 4.474 

- a meson field constants g2 = -10.431 1/fm, g3= -28.885 

The relativistic Hartree equations are solved by iterations: one starts with some 
estimate of the meson and electromagnetic fields, then solving the Dirac equation one 
finds the Dirac spinors. They give the densities and currents as sources for the Klein- 
Gordon equations for the meson fields. After their solution the new set of the meson 
and electromagnetic fields is found as the starting point for the next iteration. When 
the selfconsistency is achieved, Hartree-Bogolubov wave functions \|/p("^) of protons and 
neutrons are used for evaluation of the mean values of operators in interest. We assume 
here the product of the BCS-type functions 



I*) = Yl{u„ + v^a^at^)\vac) 

u>0 



(7) 



for protons and neutrons as the ground state wave function of nucleus. 

To get the strength of the pairing interaction for nuclei from very different regions, 
we have taken simply the experimental energy gaps (A) and the lowest in energy Z (or 
N) single-particle levels when solving the BCS equations. Such a procedure is justified 
in our calculations because the quantities which we evaluate depend rather weakly on the 
choice of the pairing force. All A-s are extracted from the experimental mass differences 



taken from the Wapstra-Audi tables ^^ but for the isotopes for which the experimental 
data do not exists we have used the estimates: Ap(„) = 12/vA MeV. 
The monopole moments of protons and neutrons distributions are 

Qf^ = (^lE^'l^)"^"^- (8) 

The mean square radii (MSR) are defined 

{r\ = f , {r\ = f . (9) 

and the root mean square radii (RMSR) are 

Tpin) = ^ir^)p(n) = \i^Rp{n) ■ (10) 

In equations (9-10) we have neglected corrections originating from the center of mass 
motion. For heavier nuclei which we discuss these corrections are small. The quadrupole 
moments of proton and neutron distributions are given by 

g^W = (v^|5]2r2p2(cos^,)|vl/)P("), (11) 

V 

where P2 is the Lagrange polynomial of the order 2. The quadrupole deformation param- 
eter of proton or neutron distribution is approximately equal to 

(12) 

We have evaluated also the binding energies of nuclei (-Brmf), the reduced electric quadru- 
pole transition probabilities {B{E2)), the proton and neutron separation energies (S'p(„)). 
The results are compared with the experimental data taken from Refs. [1, 2, 14]. 

3 Results 

Various functions of radii and density moments were investigated in order to extract the 
isospin dependence of proton and neutron density distributions. The calculations were 
performed for the stable nuclei along the (3 stability line and for all potentially existing 
isotopes with Z=50, 56, 82, 94 and isotones with A^=50, 82, 126. Schematic view in the 
{N, Z) plane of the nuclides discussed in the paper is presented in Fig. 1. For the nuclei 
out of /9-stability line the calculation was made until the proton or neutron drip line was 
reached. 

In Fig. 2 the set of results concerning the /5-stable nuclei with AO < A < 256 is 
presented as a function of mass number A. For each A value only one isotope with the 
smallest mass is chosen. 




The upper left figure shows the difference between the binding energy -Brmf of 
a nucleus calculated by the RMFT with the NL-3 parameters set and its experimental 
value -Bexp taken from [^ . Here the maximal error of the RMFT estimates is 5 MeV only 



(for A = 150). 

The upper right figure presents the reduced electric quadrupole transition probabil- 
ities B{E2) obtained by the RMFT (solid line) compared to the experimental data 0] 
(circles). The agreement is rather good in spite of too large theoretical predictions of the 
nuclei about A ~ 170 and 250. 

The separation energies of neutron Sn (left hand side) and proton Sp (right hand 
side) evaluated within the RMFT agree with the experimental data |T^ (circles) very 
well. 

The neutron r„ and proton r^ root mean square radii calculated by the RMFT 
presented in the lowest left hand side figure are also close to the experimental data (crosses 
for protons |I| , circles for neutrons |0]). The both radii are slightly different from each 
other, what is connected with the different in size and in deformation proton and neutron 
density distributions. 

The difference of the proton and neutron quadrupole deformations obtained in the 
RMFT calculation is shown in the lowest right part of Figure 2. One can see that for 
same nuclei it exceeds 0.03. 

The next seven figure sets (Figs. 3-9) show similar results as in Fig. 2 but for 
the isotope and isotone chains presented in Fig. 1. The experimental data (crosses for 
protons, circles for neutrons) are drawn, if they exist. The results should illustrate the 
goodness of the RMFT also for the nuclei outside /3 stability hne. Using these results we 
can extract the isospin dependence of proton and neutron radius constants. 

In Fig. 3 the same quantities as in Fig. 2 for Sn isotopes with A^ = 50—86 are drawn. 
All Sn nuclei are almost spherical so the quadrupole deformation of proton and neutron 
distribution is (lowest right hand side picture). Also the reduced electric quadrupole 
transition probabilities B{E2) proportional to the square of quadrupole moment are zero 
in the theory and the experiment as well (upper right hand side picture). The results 
both for proton and neutron radii are well confirmed by the experimental data what can 
be seen in the lowest left figure. The RMFT radii of neutron and proton distributions 
differ from each other up to 0.5 fm for the heaviest isotopes. The experimental separation 
energies are very well reproduced by the RMFT calculation (middle pair of figures). The 
binding energies obtained by the RMFT differ only up to 2.5 MeV from the experimental 
ones (upper left hand side picture). 

The Ba isotopes presented in Fig. 4 are slightly deformed. One can see in the lower 
right figure that the neutron and proton deformations are slightly different. The RMFT 
reproduces the experimental kink in neutron separation energy Sn around neutron magic 
number (A^ = 82) 

The Pb isotopes are presented in Fig. 5. The equilibrium deformations are largest 
for he isotopes around A = 195. The proton and neutron radii differ from each other much 
especially for the lighter isotopes. Unfortunately the RMFT estimate of the neutron root 
mean square radius for ^°^Pb is pretty far from the experimental value [^. 



The Pu isotopes, shown in Fig. 6 behave quite regularly in spite of the relative 
large deformations and B{E2). The quadrupole deformations of proton and neutron 
distributions differ by about 0.025 for all isotopes. This difference is about 10 % of the 
average deformation of these nuclei. 

In spite of the constant Z values the proton radius in Sn, Ba and Pb isotopes grows 
with the neutron number. It is not the case for the isotones shown in the next 3 figures: 
7-9, where the neutron radii stay almost constant when number of protons Z grows. 

In Fig. 7 the A^ = 50 spherical isotones (Ni-Sn) are presented. The errors in the 
binding energy are rather small. Due to the spherical shape of these nuclei no differences 
in proton and neutron deformations are observed. Contrary to the neutron radius the 
proton radius grows here rapidly with A. 

Also the A^ = 82 isotones (Cd-Hf ) shown in Fig. 8 have an almost constant neutron 
density distribution radius. This conclusion is also confirmed by the results presented in 
Fig. 9, for A^ = 126 isotones of Hg-U. 

In order to extract the average isospin dependence of the radius constant from the 
root mean square radii evaluated within the RMFT we have divided the RMSR values 
by the factor J3/5A^^^ (see Eq. 10). As one can see in Figs. 10a and 10b the radius 
constants of protons (crosses) and neutrons (circles) are gathered around 1.25 fm and 
1.20 fm respectively as functions of the relative neutron excess / = (A^ — Z)/A. In order 
to analyze better the results one has to remove the influence of deformation on nuclear 
radii. We have renormalized them to the sphere using the volume conservation rule. It 
was made approximately by dividing the nuclear radii constants by factor g: 

9{f32) = l + ^Pl (13) 



This renormalized quantities Rq = r/(w3/5 ■ g{f32)) are shown on Fig. 10a, b for 
Ba and Pb isotopes respectively. Now one can see that the renormalized proton radii 
(diamonds) decrease almost linear with I while the neutron ones (dots) increase. 

The above investigation has convinced us that the microscopic radius constant rg'" 
must depend also on the isospin /. The additional dependence on 1/A was noticed when 
comparing the results for different isotope chains. In Figs. 11, 12 the results are shifted 
by the terms proportional to 1/A in order to see better the average / dependence of the 
radii. 

In Fig. 11 we present the proton and in Fig. 12 the neutron renormalized radii 
constant for all the isotope and isotone chains discussed in this paper (a) and for the 
/5-stable nuclei (b). The results for protons are shifted by the term 0.8 /A fm while the 
neutron estimates by —3.3/ A fm. 

The dashed lines present average behavior of calculated radii constants. It was found 
by minimization of the mean square deviations. The following formulae for proton and 
neutron radii: 

R^ = 1.237 (l - 0.157 -^ 0.646 ■ -) A^^^ fm 



and (14) 

R^ = 1.176 (l + 0.250 -— + 2.806 ■ -\ A^^^ fm 

approximate the results obtained within the RMFT with the NL-3 set of the parameters. 
One can beheve that the parameters of these formulae, found for the representative nuclei 
all over the periodic system can describe well the average trend of the results obtained 
within the RMFT. 

The following formula approximates the RMFT results for the nuclear charge radii: 

i?f = 1.241 (l - 0.154 -— + 0.580 ■ -\ A^^^ fm . (15) 

This equation, obtained only by the analysis of the RMFT calculation results, has not 
much different parameters to those from equation (3) fitted to all available experimental 
data. The coefficient of (A^ — Z)/A term is only slightly smaller than in the phenomeno- 
logical formula (3). But the parameter at the term 1/A in (15) is almost 3 times smaller 
than in Eq. (3). It is mostly due to the fact that we have analyzed the nuclei with A> AO 
while the Eq. (3) was obtained by fitting the experimental data for all nuclei with A > 12. 

From Figs. 11 and 12 we can learn that the average formulae (14) work properly 
only for nuclei with A > 60. Also some nuclei with / ~ 0.2 don't fit to the average 
formulae. It is caused by non quadrupole deformations (octupole?) of these nuclei which 
were not included in the present analysis. 

As far as the ratio of proton to neutron radii (or RMSR) is concerned there is almost 
no influence of deformation on the results. It is caused by the similar or identical shapes 
of proton and neutron density distributions. We have found a short formulae for this ratio 
similar to those for radii constants (14,15). The dependence of proton to neutron radii 
ratio on / is shown in Fig. 13. We have shifted the ratio by the term —3.3/A in order 
to remove the influence of the 1/A dependence. Dashed lines in Fig. 13 represent the 

formula: 

r / N — Z 1 \ 

^ = 1.048 1 - 0.364 3.148 ■- (16) 

fitted by the least square fit to the Vp/vn RMFT ratios for all nuclides shown Fig. 1. 
In Fig. 13a the radii ratios of these chains of isotopes and isotones are presented while 
in Fig. 13b the corresponding results for /9-stable nuclei with A > 40 are plotted. The 
results for nuclei with A > 60 from are reproduced well by formula (16). 

Eq. (16) can be used to estimate the neutron RMSR from the charge one. It could 
be useful because it exists over 250 measured charge RMSR and only few experimental 
ones for the neutron RMSR. One can rearrange Eq. (16) using the relation between the 
charge and proton mean square radii 

r^^ = rl + 0.64 fm^ . (17) 

The term 0.64 fm^ originates from the final size of the proton. Finally one gets: 



d - 0.64)V2 



1.048 (1 - 0.364 ■ {N - Z)/A - 3.148/A) 



f m , 



In Fig. 14 all known experimental RMSR (errorbars) for neutrons are compared 
with the RMFT predictions (crosses) and with the results estimated from experimental 
charge RMSR 0] using Eq. (18) (circles). All results agree very well with each other, so 
we hope that the formula (18) can be used by the experimentalists to foresee the neutron 
distribution radii for other nuclei as well. 

4 Conclusions 

The following conclusions can be drawn from our calculation 

1. The renormalized to sphere RMFT proton and neutron distributions radii i?o 
depend almost linear on neutron excess / = (A^ — Z^j A. Rq decrease with /, while 
Rq increase what was suggested in |^ and obtained in by the analysis of the 
results obtained within the HFB calculation with the Skyrme forces. 

2. The parameters of the formulae for the charge radii are similar as in the phenomeno- 
logical formula in [§. However the dependence of the RMFT radii on / is about 
20% weaker than in the phenomenological formula which describes the global exper- 
imental trend. It means that the parameters NL-3 of the RMFT should be slightly 
changed. 

3. The term ~ A~^ in the formula for Rq is needed in order to reproduce the average 
MSR values obtained within the RMFT. Its role is especially important for lighter 
nuclei. 

4. The ratio of the proton radius to the neutron radius is a smooth function of I and 
A~^ and it could be very well described by the simple formula (16). Using this 
global dependence we have written the phenomenological formula (18) which allows 
to foresee the magnitude of the neutron radius when the experimental charge radius 
is known. The prediction power of the formula (18) is not worse than that of very 
advanced microscopical calculations based on the RMFT or the HFB-Skyrme model 

i- 

The obtained formulae for the Rq and Rq radii will be used to develop the liquid drop 
like model, which will depend on the different proton and neutron density distributions. 



i.e. different radii and deformations. Such investigations are in progress now ||T5|. 
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Figure captions 

1. Schematic view in the {N,Z) plane of the nuchdes discussed in the paper. For 
the nuclei out of the /^-stability line the calculation was made until the proton or 
neutron drip line was reached. 

2. The results obtained within the RMFT+BCS model with the NL-3 parameters [|r^ 
for the /3-stable nuclei compared with the experimental data as a functions of mass 
number A. In the diagrams are presented: the difference between the calculated 
-Brmf and experimental -Bgxp binding energy (upper-left figure), reduced electric 
quadrupole transition probabilities B{E2) (upper-right figure), neutron and proton 
separation energies Sn, Sp (middle figures), neutron and proton root mean square 
radii rn,rp (lower-left figure), differences between proton and neutron quadrupole 
deformation parameters /Jf — P2 (lower-right figure). 

3. Same as in Fig. 2 but for the Sn isotopes {Z = 50). 

4. Same as in Fig. 2 but for the Ba isotopes {Z = 56). 

5. Same as in Fig. 2 but for the Pb isotopes {Z = 82). 

6. Same as in Fig. 2 but for the Pu isotopes (A^ = 94). 

7. Same as in Fig. 2 but for the isotones with A^ = 50. 

8. Same as in Fig. 2 but for the isotones with A^ = 82. 

9. Same as in Fig. 2 but for the isotones with A^ = 126. 

10. The proton (crosses) and neutron (open circles) nuclear radii divided by A^^^ as 
function of the relative neutron excess / for Ba (a) and Pb (b) isotopes. Nuclear 
radii renormalized to sphere are plotted with diamonds and dots respectively. 

11. Renormalized to sphere proton radii constants for discussed chains of isotopes and 
isotones (a) and for all (3 stable nuclei with A > AO (b) as function of the relative 
neutron excess / shifted down by the —3.3/A fm term. Dashed line represents 
average behavior of calculated values for the chains (14). 

12. Same as in Fig. 11 but for neutron radii (shifted up by the 0.8/ A fm term). 

13. The ratios of the proton to neutron root mean square radii for discussed chains of 
isotopes and isotones (a) and for the /3 stable nuclei with A > 40 (b) as functions 
of relative neutron excess I. The ratios are shifted up by 3.3/A fm term. Dashed 
lines represent the average behavior of calculated values (16). 

14. Comparison of the experimental neutron root mean square radii [Q (errorbars) with 
the RMFT results (crosses) and the values obtained from experimental charge radii 
[0] through the new formula (18) (circles). 
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